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Abstract 

The statistics of the hght emitted by two-level lasers is evaluated on the basis of 
generalized rate equations. According to that approach, all fluctuations are inter- 
preted as being caused by the jumps that occur in active and detecting atoms. The 
intra-cavity Fano factor and the photo-current spectral density are obtained analyt- 
ically for Poissonian and quiet pumps. The algebra is simple and the formulas hold 
for small as well as large pumping rates. Lasers exhibit excess noise at low pumping 
levels. 
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1 Introduction 



We consider tlie liglit statistics of two-level atom lasers. Analytical expres- 
sions for the photo-detection spectral density, ^(fi), and the Fano factor, 
^ = var(m)/ (m), where m denotes the number of photons in the laser cav- 
ity, var(m) the variance of m, and (m) the average value of m, are derived. 
Numerous paper have addressed this problem before, e.g. 0, |], H, |, ^, ^ ||. 
However, explicit formulas that are not restricted to high-pumping levels and 
do not neglect spontaneous decay do not seem to be available. The concepts 
on which rests our generalized rate-equations approach have been explained in 
detail in tutorial papers |^, |T0| and applied to three- and four-level lasers |TT|] . 



This method enables one to obtain analytical expressions straightforwardly, 
which coincide with Quantum-Optics results whenever a comparison can be 
made, even in the case of sub-Poissonian light statistics. The novelty of the pa- 
per resides in the fact that a simple formula is found that holds irrespectively 
of the pumping rates. Similar formulas were previously obtained for four-level 
lasers []TT| , but they were too involved to be written down. The present paper 



may be viewed as an illustration of that previous work, with formulas that are 
easier to appreciate because of their simplicity. 

Rate equations treat the number of photons in the cavity as a classical random 
function of time. The light field is quantized as a result of matter quantization 
and conservation of energy, but not directly. Rate equations should be dis- 
tinguished from semi-classical theories in which the optical field is driven by 
atomic dipole expectation values. The present theory rests instead on the con- 
sideration of transition probabilities as in []T2|. Every absorption event, even 



if it occurs in the detector, reacts on the number of photons in the optical 
cavity. Our method allows the use of Monte-Carlo simulations [jll, 13 . 



The active medium considered is a collection of identical atoms. Transitions 
are supposed to be strongly homogeneously broadened so that atomic polar- 
izations may be adiabatically eliminated. The two-level scheme is obtained as 
a limit case of the three-level one (see Fig. 1), with the levels labeled |0), 
and 1 2) in increasing order of energy. Level separations are supposed to be 
large compared with k-^T, where T denotes the atomic temperature and 
the Boltzmann constant, so that thermally-induced transitions are negligible. 
The classical 2-level scheme is a valid one in the limit in which the lower decay 
time 1/pd tends to zero. In that limit, the |1) population tends to zero and 
stimulated absorption from |1) to |2) does not occur. For concreteness, we may 
think of the pumping as resulting from the injection of atoms in the excited 



state inside the cavity, as in Scully and Lamb |TJ]. If the atoms are injected 



regularly the fluctuations of the pumping rate vanish = in our notation). 
The probability per unit time that an electronic transition from level |2) to 
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Fig. 1. Level scheme of a three-level laser. The two-level laser obtains when 1/pd — > 0. 

|1) occurs is taken as equal to m+ 1, where m denotes the number of photons 
in the cavity. This amounts to selecting a time unit for which the laser gain 
(7 = 1. Spontaneous decay from level |2) to level |1) is allowed with probability 
7. This decay may be either non-radiative or involve radiation into other elec- 
tromagnetic modes, besides the one of interest. Photons are absorbed with 
probability am, where a denotes a constant, the absorbing atoms residing 
most of the time in their ground state. These absorbing atoms model the 
transmission of light through mirrors with subsequent absorption by a detec- 
tor. Provided detection is linear and refiectionless, it is immaterial whether 
absorption occurs inside or outside the optical cavity. For simplicity, internal 
absorption is neglected. 



Golubev and Sokolov first pointed out in |1984| that if the laser pumping 
rate does not fluctuate and detection is ideal, the detected current should 
not fluctuate at low Fourier (or baseband) frequencies. This observation was 
made independently by Yamamoto et al. [|I|. These latter authors furthermore 
proposed |1986| a practical way of realizing quiet electrical pumps, namely 
a battery in series with a cold resistance, much greater than the laser diode 
dynamic resistance [|16i. This scheme was implemented for the first time in 



19871 . The pumping rate is supposed in the present paper to be fully prescribed 



through its average value J and the spectral density aj = of its fluctuations. 
It would be straightforward, however, to account for a dependence of the 
pumping rate on the atomic populations, as is required in the case of optical 
pumps ITT 



We first consider the steady-state regime and subsequently derive formulas for 
the photo-detection spectrum and the intra-cavity field Fano factor. 



2 Steady-state 



Let rij, j — 0, 1, 2, denote the number of atoms in state j. In the hmit 1/pd — 0, 
we have rii — 0, and thus 

N = no + n2. (1) 



Let J' denote the pumping rate, TZ the net stimulated rate, S the spontaneous 
decay rate from the upper to the lower working levels, and Q the photon 
absorption rate. The steady-state conditions read 

j^n + s, Q = n, (2) 

where 

J — J, TZ — (m + l)n2, S — 7^2, Q — am. (3) 



m 



Equations (1), (2) and (3) provide the steady-state atomic populations rij and 
photon number m as 

a/(J + 7+1) -4J7-7-1, 712 = ^—, no = A^ — , 4 

y ^ ' m-|-l m-|-l 

where J = J/a. 

When 7 = 0, Eq. (4) leads to 

- aJ aJ 

m — J, 712 = — , nn — JS — . (5) 

J+1 J+1 



3 Photo-detection spectrum 



Within the weak-noise approximation atomic and photonic populations split 
into steady-state values and fluctuations. For example the instantaneous pho- 
ton number m may be written as m = (m) + Am, where (m) needs not be not 
distinguished from the steady-state value of m. Similarly, rates get split into 
steady-state values and fluctuations consisting of a deterministic function of 
the population fluctuations and uncorrelated Langevin "forces" . For example 
J splits into J = {J') and the sum, AJ, of a deterministic function of pop- 
ulations and a Langevin force j{t) expressing the jump process randomness. 
Thus 

J = J + AJ, n = R + AR, S = S + AS, Q = Q + AQ. (6) 
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where a first-order variation of tlie expressions in Eq. (3) gives 



A J = J, (7a) 

Ai? = (m + 1) An2 + n2Am + r, (7b) 

AS = 7An2 + s, (7c) 

AQ = aAm + q. (7d) 

Since tfie number of atoms is supposed to be constant we liave 

Ano + An2 = 0. (8) 

At some Fourier frequency fl, [3] the generahzed rate equations read ^ 

iQAm = AR- AQ, (9a) 

inAno = AR + AS- A J, (9b) 

ifiAna = A J - Ai? - AS. (9c) 



Equations (6)-(9) are solved for AQ, yielding a linear combination of the 
Langevin forces 

AQ= J2 c.;^, (10) 

z&{j,q,r,s} 

where the are complex and frequency dependent coefficients. The normal- 
ized photo-current spectral density is thus 

y{n) = — Yl \czf(T,, (11) 

am ^ f 1 

where denotes the spectral density value of the Langevin noise source z, 
equal to average rates 

aj = ^J, CTg = am, ar = {m + 1) n2, CTg = 'yn2. (12) 

^ accounts for a poissonian pump if ^ = 1, while a perfect quiet pumping 
is obtained at ^ = 0. When these expressions are introduced in Eq. (11) an 
analytical expression of is obtained 



^ The Fourier angular frequency is called here 'frequency' for short. The frequency, 
as it is usually defined, is obtained by dividing by 2tt. 
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with 



^ (VC) = {n2{3m + 1) — am) + n2(3m + l)7m — amm^ + J{m + 1)^.^, 

(14a) 

= Vl'^ + ((m + a - + 2 (7^2 - am)) fi^ + (7^2 - amf . (14b) 

where m stands for m + 7 + 1. S^[VL) — 1 is the ratio of a polynomium of 
degree 1 in in the numerator and of degree 2 in VL^ in the denominator. A 
few observations can be made from Eqs. (13)-(14) 

• Well above threshold and with 7 = 0, implying that m S> 1 and J ^ 1, 
Eqs. (13)-(14) simplify to y = [a^^ + Vt^)/{a'^ + Vt^)^ as is well-known, see 

e.g. i; 

• lim^^oo -$^(^) = 1. That is, the light statistic is Poissonian at high Fourier 
frequencies irrespective of the pumping and laser parameters; 

• ^(0) = ^ for an ideal laser (7 = 0), i.e., laser light spectrum at zero- 
frequency corresponds to the value of the pump ^-parameter; 

• For usual laser parameters there is a value of VL that maximizes 
according to Eq. (13). 

The features of y{Q) just discussed are illustrated in Fig. 2. We have chosen 
two- level laser parameters corresponding to a thresholdless laser (7 = 0), 
N = 10^ atoms, and cavity losses a = 6.32 identical to those considered for 
three- and four-level lasers in [|TT|]. Figs. 2a and 2b verify that two- level lasers 



are noisy even when spontaneous decay is negligible. This occurs at a Fourier 
frequency that depends on the pumping level. The maximum spectral-density 
value is a function of the number m of photons stored in the cavity. The curves 
given in Figs. 2c and 2d demonstrate that the laser is most noisy when m is 
on the order of unity. At small m-values there is no phase-coherence between 
successive wave trains emitted because they do not overlap, and the light is 
thermal-like. At high m-values a regulation mechanism takes place. 



4 Fano factor 



The intra-cavity photon statistics is characterized by the Fano factor ^ = 
vai{m)/avem. The variance of m is the integral over frequency of ^Am(^); 
the normalized spectral density of Am. Similar to Sec. 3, ^Am(^) is obtained 
by solving Eqs. (6)-(9) for Am instead of AQ 
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(c) (d) 

Fig. 2. Spectral density of two-level lasers with 7 = 0, = 10^ and a = 6.32. 
Upper curves are 3D plots of =5^(0, J). Lower curves give =5'^ax ^ function of ui, 
the photon number within the cavity. At a fixed J, ^max is the maximum value of 
=5^(0, J) when Q is varied: (a) and (c) ^ = 1, Poissonian pump, (b) and (d) = 0, 
quiet pump. 

where 

y{n) = {n2{m + 1) + am) + n2{m + l)7m + amm^ + J{m + if^, (16) 
and with ^{Q) still given by Eq. (14b). 
The Fano factor is thus 

/°° HO 
-00 ztt 

for which a lengthy analytical expression has been obtained using symbolic 
computations. Well above threshold and with 7 = 0, this expression yields 
^ = (1 + ^) /2, in agreement with [Q. 



7 




Fig. 3. Fano factor of two-level lasers with parameters N = 10^ and a = 6.32 as 
a function of the pump rate J and the spontaneous emission rate 7. (a) ^ = 1, 
poissonian pump, (b) ^ = 0, quiet pump. (I) 7 = 0, (II) 7 = 6.32, (III) 7 = 63.2, 
(IV) 7 = 632, (V) 7 = 6325. 

The Fano factor of two-level lasers is represented as a function of the pump- 
ing rate and the spontaneous decay rate in Fig. 3. When 7 = a bump is 
observed on at J ~ 10, related to the bump observed on ^ in Fig. 2: light 
statistics is nearly that of thermal light. The bump sharpens as 7 increases 
with a maximum occurring at threshold. Below threshold the light statistics is 
thermal-like, so that ^ ~ (m) + 1. Above threshold Poissonian light statistics 
(^ ^ 1) is being approached, for the case of a Poissonian pump. 



5 Conclusion 



Two-level lasers have been considered within the framework of generalized 
rate-equations. Fully analytical expressions of the photo-current spectral den- 
sity and intracavity field Fano factor were obtained. It was shown that thresh- 
oldless two-level lasers become noisy when the average number of photons in 
the cavity is on the order of unity. Sub- Poissonian light statistics is accurately 
described. The present paper should be considered as being partly tutorial, in 
that it provides more accessible formulas than does a previous paper dealing 



with four- level atom lasers 11 



References 

[1] Y. Yamamoto, S. Machida, and O. Nilsson, Phys. Rev. A 34, 4025 (1986). 
[2] T. A. B. Kennedy and D. F. Walls, Phys. Rev. A 40, 6366 (1989). 
[3] A. M. Khazanov, G. A. Koganov, and E. P. Gordov, Phys. Rev. A 42, 
3065 (1990). 



8 



[4] H. Ritsch, P. Zoller, C. W. Gardiner, and D. F. Walls, Phys. Rev. A 44, 
3361 (1991). 

[5] T. C. Ralph and C. M. Savage, Phys. Rev. A 44, 7809 (1991). 
[6] M. I. Kolobov, L. Davidovich, E. Giacobino, and C. Fabre, Phys. Rev. A 
47, 1431 (1993). 

[7] R. B. Levien, M. J. Collett, and D. F. Walls, Phys. Rev. A 47, 5030 
(1993). 

[8] G. A. Koganov and R. Shuker, Phys. Rev. A 63, 015802 (2000). 
[9] J. Arnaud, Opt. Quantum Electron. 27, 63 (1995). 
[10] J. Arnaud, Opt. Quantum Electron. 34, 393 (2002). 
[11] L. Chusseau, J. Arnaud, and F. Philippe (2002), [quant-ph/0203029 . 
[12] R. Loudon, The Quantum Theory of Light (Oxford University Press, Ox- 
ford, 1983). 

[13] L. Chusseau and J. Arnaud, Opt. Quantum Electron. (2002), to appear, 

|quant-ph/010507£ . 
[14] M. O. Scully and W. E. Lamb, Jr., Phys. Rev. 159, 208 (1967). 
[15] Y. M. Golubev and I. V. Sokolov, Sov. Phys.-JETP 60, 234 (1984). 
[16] S. Machida, Y. Yamamoto, and Y. Itaya, Phys. Rev. Lett. 58, 1000 

(1987). 



9 



